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FOREWORD 


This report covers work completed during the period February 1, 1978 
through July 31, 1978 on the research project "Influence of Precursor Heating 
on Nonequilibrium Viscous Flow Around a Jovian Entry Body." The work was 
supported by the NASA/Langley Research Center (Aerothermodynamics Branch of 
the Space Systems Division) through research grant NSG 1492. The grant was 
monitored by Dr. Randolph A. Graves, Jr. of SSD-Aerothermodynamics Branch 
(Mail Stop 366) . 
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EFFECTS OF PRECURSOR HEATING ON CHEMICAL AND RADIATIVE NONEQUILIBRIUM 
VISCOUS FLOW AROUND A JOVIAN ENTRY BODY 

By 

S. N. Tiwari 1 and K. Y. Szema 2 
SUMMARY 

The influence of precursor heating on viscous chemical nonequilibrium 
radiating flow around a Jovian entry body is investigated. Results obtained 
for a 45-degree hyperboloid blunt body entering Jupiter's nominal atmosphere 
at zero angle of attack indicate that the nonequilibrium radiative heating 
rate is significantly higher than the corresponding equilibrium heating. 

The precursor heating, in general, increases the radiative and convective 
heating to the body, and this increase is slightly higher for the nonequilibrium 

I 

conditions . f 

'i 

•i 

NOMENCLATURE 

C^ mass fraction of species i in the shock layer, p^/p 

J^ mass diffusion flux of species i 

k Boltzmann constant 

k^,r backward rate constant 

k c ,r equilibrium rate constant 

k^,r forward rate constant 

Le Lewis number 

M* molecular weight of mixture 

N number of moles 

1 Professor, Department of Mechanical Engineering and Mechanics, Old Dominion 
University, Norfolk, Virginia 23508. 

2 Graduate Research Assistant, Old Dominion University Research Foundation, 

P. 0. Box 6369, Norfolk, Virginia 23508. 




*3 


n 




P 

Pr 


*^R 

r 


r 

s 

R* 


R* 

s 

s 

T 


T* 

ref 


u 

u' 

v 

V* 

X 

y 

a 

£ 

e 

n 

K 



coordinate normal to the bow shock, n*/R* 

pressure, p*/(p* V* 2 ) 

Prandtl number, y* C*/k* 

P 

net radiant heat flux, q£/(p* V* 3 ) 

radius measured from axis of symmetry to a point on the body 
surface, x*/R* 

radius measured from axis of symmetry to a point on the bow 
shock, r|/R* 

universal gas constant . 

radius of the body 

body nose radius (same as R*) 

T1 

radius of the bow shock 

coordinate along the bow shock, s*/R^ 

) 

temperature, T*/T* ef | 

reference temperature, V^/C^ | 

velocity tangent to body surface, u*/V* 
velocity tangent to bow shock, cm/sec 

velocity normal to body surface, u*/V* i 

velocity normal to bow shock, cm/sec 

\ 

coordinate along the body surface, x*/R^ 
coordinate normal to the body surface, v*/Rj$ 
shock angle defined in figure 1 
Reynolds number parameter or surface emittance 
body angle defined in figure 1 
transformed y coordinate, y/v 

5 

body curvature, tc*R^ 
viscosity of mixture, U*/u£ e £ 
reference viscosity, u*(TJ e £) 
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coordinate along the body surface, £=x 


p density of mixture, p*/p* 

00 


Subscript 

i 

s 

w 

09 

V 


ith species 
shock value 
wall value 

free* stream condition 
radiation frequency 


1. INTRODUCTION 

In two previous studies (refs. 1 and 2), the authors have investigated 
the influence of precursor heating on the flow phenomena around a Jovian entry 
body. In these studies, the precursor as well as shock-layer gas was assumed 
to be in chemical nonequilibrium. While for most entry conditions envisioned 
for Jupiter's atmosphere it may be reasonable to assume chemical equilibrium 
in the precursor region, this assumption could not be justified logically 
for the shock layer. This is because (aside from inclusion of the precursor 
effect) there exists strong evidence of a significantly large chemical 
nonequilibrium layer (closer to the shock) in the shock layer for Jovian 
entry conditions (refs. 3 to 6). Also, sutdies on heating of high-speed 
vehicles entering other planetary atmospheres (refs. 7 to 10) have indicated 
that, under certain conditions, the radiation from the nonequilibrium region 
of the shock layer was considerably higher than that predicted by an 
equilibrium analysis. 

The purpose of this study, therefore, is to investigate the influence of 
precursor heating on the flow phenomena around a Jovian entry body under the 
chemical nonequilibrium in the shock layer. The reaction scheme (which 
includes hydrogen dissociation, electronic excitation, and ionization by 
atom-atom and atom-electron collisions) proposed by Leibowitz (ref. 3) is 
used in this study. Ionization of hydrogen is initiated by atom-atom collisions 
which produce atoms in electronically excited states. These are rapidly 
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ionized by additional collisions. A simplified two-step model is assumed in 
which hydrogen in the first excited state is produced by collisions with 
hydrogen and helium atoms and is then immediately ionized by a subsequent 
collision. The hydrogen is then ionized by collisions with electrons. 

Molecular hydrogen is first dissociated by collisions with itself and helium. 
The hydrogen atoms produced in the reaction also participate in the dissocia- 
tion reaction as do ions and electrons as they are formed. The solutions of 
the governing viscous shock layer and inviscid precursor equations are 
obtained by employing the numerical procedures outlined in references 1 and 2. 

2. ANALYSIS 

The physical model and coordinate system for a Jovian entry body is 
shown in figure 1. As discussed in references 1 and 2, the entire flow field 
ahead of the body is divided into three regions: the free stream, the precursor 
region (preheating zone), and the shock layer. The basic governing equations 
and boundary conditions for the present problem are essentially the same as 
those given in reference 2. There are, however, a few definite changes in the 
governing equations because of the chemical nonequilibrium condition, and 
these are discussed briefly in this section. 

2.1. Governing Equations for the Precursor Region 

The precursor region is considered to be thin, and the flow in this 
region is assumed to be steady, inviscid, and in chemical equilibrium. The 
governing equations for this region, therefore, are the same as those given 
in reference 2. Because of the chemical nonequilibrium in the shock layer, 
however, the conditions at the outer edge of the shock are found to be 
different from those obtained in reference 2. It should be pointed out 
here that, for the coupled precursor- shock layer flow phenomena, the 
conditions at the outer edge of the shock are obtained through iterative 
procedures . 


2.2. Shock-Layer Equations 

The flow in the shock layer is assumed to be axisymmetric, steady, 
laminar, viscous and compressible. The shock layer gas is assumed to be in 
local thermodynamic equilibrium but in chemical nonequilibrium. The basic 
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governing equations (continuity, x -momentum, y-momentum, and energy), there- 
fore, are the same ones given in reference 2. The species continuity equation, 
however, is different, and this is given by the relation 

P [(u/r) (3C./3x) + vCaCj/ay)] = • (e 2 /Tc) 

*[o/3y) crej.)] ( 1 ) 

where 

T » 1 + yic, ? = r + y cos 6, e * uj ef /(p* V* RJ) l/Z , 

J. = (u/Pr) Le (3C ./3y) 

and w i represents the rate of production of chemical species in the shock 
layer. The equation of state proposed by Zoby et al. (ref. 11) is valid 
only for chemical equilibrium conditions. In the case of chemical nonequi- 
librium, the equation of state is given by the relation (ref. 12) 

P*v* - (I N.)R*T* (2 ) 

where is the number of moles for the ith species. This result is 
reminiscent of the thermal equation of state for a perfect gas. The sum in 
the parentheses, however, is not a constant since the total number of moles 
changes as th > chemical balance changes. Other shock-layer equations are 
the same as those presented in reference 2. 

As in reference 2, no-slip and no temperature- jump conditions are used 
at the body surface (wall). Thus, % * v w * °» and the wall temperature is 
either specified or calculated. Furthermore, it is assumed that the chemical 
species are in equilibrium at the body surface. The conditions in front of 
the shock are obtained from the precursor solutions; the conditions immediately 
behind the shock are obtained by using the Rankine-Hugoniot relations. Two 
assumptions can be made about the molecular hydrogen entering the 3 hock layer 
immediately behind the shock. One criterion is to assume that chemical 
reactions are "completely frozen" and that initial composition of hydrogen 
just behind the shock corresponds to the free-stream value. The second 
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criterion is to consider that all hydrogen moledules have been dissociated 
immediately behind the shock. This is referred to as the "half frozen" 
condition. Both conditions are used in this study. 


2.3. Rate of Species Production in the Shock Layer 


The reaction scheme describing important collisional processes in 
hydrogen-helium ionizing shock waves has been modeled by Leibowitz 
(refs. 3 and 6) after the results of argon ionization studies. Eleven 
separate reaction steps describe the dissociation of molecular hydrogen, 
excitation of electronic states of hydrogen and helium, and ionization of 
the atomic hydrogen and helium by collisions with atoms and electrons. A 
complete discussion on these reaction schemes is available in reference 3, 
and a brief discussion is presented here. The eleven reactions and corresponding 
rate constants are given in table 1. 


In a complex gas mixture containing a total of l species of x., 
capable of undergoing m elementary chemical reactions, the chemical 
equation for the general elementary reaction r can be written as 
(ref. 12) 



(3) 


where a^ r and b^ r are the stochiometric coefficients appearing on the 
left and right in the reaction r. By applying the principle of detailed 
balancing, the backward rate constant, k^ , is obtained by dividing the 
forward rate constant, k- , by the equilibrium constant k . The 

CjT 

equilibrium constants for hydrogen, helium species were obtained by Leibowitz 
(ref. 3) by using the thermochemistry program of Horton and Menard (ref. 13) . 
These constants were obtained also in a separate study by Zoby et al. 

(ref. 14) and are expressed here as 

k^ - 4.05 x 10" 9 T 3 / 2 exp (-1.578 * 10 5 /T) (4a) 

k H + - 1.62 x io ~ 8 T 3 / 2 exp (-2. 853 x io 5 /T) (4b) 

k H - 3.7[l - exp(-1.5 x 10®/T 2 )] exp(-S2340/T) (4c) 
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The total rate of change in x^ is given by the relation (ref. 12) 



m 


L Cb i,r- i i,r’ k f,r 
r=l 


fj (Xj) a i,r 
i-1 


k 1 II (x.) b i,r 

c ' r i-1 1 


( 5 ) 


This is the general rate equation for a complex gas mixture . The rate of 
production of chemical species, w^, now can be expressed by 

w. - M* (dx*/dt) (RJ/p* v*) (6) 

Equations (1) to (6), along with other fluid mechanical equations and the 
electron energy equation, are solved numerically to obtain the concentration 
of all species. In order to have a reasonable rate of convergence in the 
numerical scheme, however, it is important to express the rate of production 
term in a proper form. This is accomplished by splitting w^ into two 
separate contributions as (refs. 15-18) 

Wj/p - (w^ 0 - (Wj) 1 (7) 

The reasons for doing this are explained in the cited references. 

2.4. Electron Temperature 

Because of a large ratio of atom (or ion) mass to electron mass, electrons 
transfer energy rapidly by collisions with other electrons but only slowly by 
elastic collisions with atoms or ions. Consequently, a different temperature 
is given to atoms (heavy particles) and electrons in the same jr s. The 
electron temperature is obtained from the solution of the electron energy 
equation. A detailed discussion of the electron energy equation is given 
by Appleton and Bray (ref. 19). For a one-dimensional, steady shock wave in a 
H 2 -He mixture, the resulting equation can be expressed as (ref. 3) : 


3[*]toe “ Te) - * ^lr + ^3 * ^3r) - 

+ R4 - R^j) - 0 


( 8 ) 
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where 


v ek ~ n k V e %k 


.V = (8k Te/nm ft ) 


Q e k * \^ ex P(-\ E ) - exp(-b k E)] (9c) 

E * (1/2)V 2 (9d) 

V = 4(k Te/2nji) l/2 (9e) 

u » n> a n^/O^ + n^) (9f) 

v = S v ek / m k (9g) 

k 

In equation (8), [e] represents the concentration of electrons, 0^ is the 
ionization energy per mole of species k, and fL and are the forward 

and backward production rates for electrons respectively. In equation (9a), 
represents the collisional frequency, n^ is the number density of 
species k, V g is the average electron velocity, and represents the 

elastic collision cross section for species k. In equations (9d) and (9e), 

E represents the relative kinetic energy, V is the relative speed, and 
u is the reduced mass. The values of coefficients A^» a^» and b^ 

appearing in equation (9c) are available in reference 3. By substituting 
equations (4a) to (4c) and (9a) to (9g) into equation (8), an implicit 
expression for the electron temperature , in terms of the species concentration 
and heavy particle temperature, can be obtained as 

Te * T - C(kx + k 3 )(x 2 - x 4 /K CjH ) * (k 2 + 1<4)(*3 - x 5/K CHe )]/*i 


where kj, k 2 , k 3 , k 4 are rate constants in table 1 and 


xi * l/(3[e]nL v R) 


(Ha) 


x 2 * [e] 


(lib) 
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W He H e ] 


(11c) 


x 3 * 

tie 

x 4 » 0 H [H + ][e] 
x 5 * 9 He [He][e] 


(lid) 

(lie) 


Use of the electron temperature is made in evaluating the radiative 
flux in the shock layer. 


2.5. Radiation Transport Model 

The radiation models used in this study are discussed in detail in 
references 1, 2, and 14. The radiative flux, q R , is calculated with the 
assumption of nonscattering medium, the tangent slab approximation for 
radiative transfer, and nonreflecting bounding surfaces. The spectral 
absorption model for the hydrogen species in the precursor region is a 
three-step model and is described in reference 1. In the shock layer, the 
absorption by the helium species is usually neglected for most Jovian 
entry conditions. In this study, a 58-step absorption model for shock 
layer gases sug'gested by Sutton (ref. 14) is used. The expression for 
the net radiative flux in the shock layer usually is given by (ref. 2): 

’r ■ % - % (l2) 

where q R r ■presents the energy transfer towards the shock and q R the 
energy transfer towards the body. For the nonequilibrium case, the radiative 
heating to the body is calculated by using the nonequilibrium species 
concentration and an effective electronic temperature. 

2.6. Thermodynamic and Transport Properties 

Thermodynamic and transport properties are required for each species 
considered in the different flow regimes. Complete information for obtaining 
these properties is given in references 2 and 18. 
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3. METHOD OF SOLUTION 


An iterative procedure has been used to couple the precursor and shock- 
layer solutions. The viscous shock-layer solutions are obtained by following 
the basic procedures suggested in references 15 to 18. The solution procedures 
for the precursor as well as shock- layer regions are described in reference 2. 

4. RESULTS AND DISCUSSION 

Information on Jupiter's atmospheric conditions is available in 
references 20 to 22. For different altitudes of entry, the free-stream 
conditions used in this study are given in references 1 and 2. The temperature 
of the atmosphere (i.e., T ) is taken to be 145 K. In the past, the nominal 
composition of the atmosphere was assumed to be 85 percent hydrogen and 15 
percent helium by mole fraction. Recently, this has been changed to 89 percent 
hydrogen and 11 percent helium (ref. 22). Illustrative results have been 
obtained by using both of these compositions. 

The entry body considered is a 45-degree hyperboloid blunt body which 
enters the Jovian atmosphere at a zero angle of attack. The body surface 
is assumed to be gray and have a surface emittance of 0.8. Unless specified 
otherwise, the surface temperature is taken to be uniform at 4,564 K. For 
the case of chemical equilibrium in the shock layer, all results were obtained 
by considering a body nose radius of R* ■ 23 cm. For chemical nonequilibrium 
conditions, however, three different nose radii (12, 23, and 45 cm) were 
considered. 

To illustrate the important features of the nonequilibrium analysis, 
most results were obtained for entry conditions which closely correspond to 
the peak heating conditions (i.e., for conditions at Z ■ 116 km). However, 
a few illustrative results have also been obtained for other entry conditions. 
Equilibrium and nonequilibrium results are presented first for variation of 
different properties in the shock layer. Results are then presented to 
illustrate the influence of precursor heating. Finally, results are 
presented for variation of different heat fluxes in the shock layer under 
the influence of both the nonequilibrium conditions and the precursor 
heating. 
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By considering initial composition of hydrogen species, just behind the 
shock, corresponding to both the "completely frozen" and "half frozen" 
conditions, nonequilibrium results were obtained for entry conditions at 
Z * 116 km and for 85 percent hydrogen nominal atmosphere. The results 
are illustrated in figures 2 to 4 as a function of the normal coordinate at 
the stagnation point. Figure 2 shows the mole concentration of different 
species across the shock layer. It is evident from this figure that molecular 
hydrogen is completely dissociated within about four percent of the total 
shock standoff distance from the shock wave. This is referred to as the 
dissociation zone (or the dissociated region). The variation in nondimen- 
sional v-velocity component and density is illustrated in figure 3. Since 
molecular weights change rapidly in the dissociated region, there is an 
increase in velocity and a decrease in density near the outer edge of the 
dissociation zone. The temperature distribution is shown in figure 4. 

It is noted that the temperature just behind the shock wave reaches a value 
of approximately 45,000 K in the completely frozen condition. After a short 
interval, however, all hydrogen molecules are dissociated and temperature 
drops to about 25,000 K. Next, ionization occurs and, as a result of this, 
temperature continues to decrease until it reaches the equilibrium value. 

From the results presented in figures 2 to 4, it is concluded that the 
half frozen and completely frozen assumptions are quite close except in 
the dissociated region near the shock wave, and that the half frozen flow 
computation is a reasonably good assumption for conditions of chemical 
nonequilibrium at altitudes near the peak heating region. Thus, all 
other results presented in this section have been obtained by considering 
only the half frozen conditions behind the shock. 

As discussed in reference 2, the shock standoff distance (for a given 
body nose radius) varies with the altitude of entry and entry velocity. 

It should be pointed out here that, in general, the shock standoff 
distance increases as the body nose radius increases. For entry 
conditions at Z ■ 116 km, equilibrium and nonequilibrium results for the 
shock standoff distance are illustrated in figure 5 as a function of the 
coordinate along the body surface. It is noted that the shock standoff 
distances for equilibrium and with radiation are considerably lower than 
for nonequilibrium and with no radiation. This, however, would be expected 
because shock-layer densities are greater for radiation and equilibrium 
conditions than for no radiation and nonequilibrium conditions. 
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Variations in chemical species across the shock layer are shown in figures 
6 to 9 for different conditions. For entry conditions at Z * 116 km, results 
presented in figures 6 and 7 show that the nonequilibrium layer is about 
25 percent of the total shock-layer thickness for the no radiation case and 
about 50 percent for the case with radiation. This is because inclusion of 
radiation results in a different temperature distribution in the shock layer. 
This point will be discussed further while presenting results for the tempera- 
ture variation. Near the wall, the mass fractions of atomic hydrogen and 
electrons are higher for nonequi librium conditions with radiation. This is 
because cold gases near the wall absorb relatively more radiative heat flux 
in the nonequilibrium case. For the case of no radiation, a comparison of 
results presented in figures 6 and 8 reveals that the nonequilibrium layer 
increases from 25 percent at Z * 116 km to about 40 percent at Z * 143 km. 

This is because density is lower at higher altitudes and, therefore, it 
will take a relatively longer time to reach an equilibrium condition. For 
Z * 116 km entry conditions, figure 9 shows the species concentrations 
for 3 different body-nose radii (12, 23, and 45 cm). These results 
indicate that the thickness (or range) of the nonequilibrium layer decreases 
with increasing nose radius. In particular, it is seen that the thickness 
is about 40 percent for R* * 12 cm, but it is only 10 percent for R* = 43 cm. 
This is because the shock standoff distance is proportional to the body nose 
radius and the relaxation time for chemical reactions is about the same for 
all cases. 

For entry conditions at Z « 116 km, absorption coefficients of the 
shock-layer gas were calculated under equilibrium and nonequilibrium 
conditions. Results obtained for conditions at a location in the shock 
layer, which is at 10 percent of the shock standoff distance from the shock 
wave, are illustrated in figure 10 as a function of electron-volts. The 
figure exhibits the trend that is quite typical of results for other 
conditions. It is seen that the nonequilibrium absorption is considerably 
higher than the equilibrium absorption for most parts of the spectral range. 

As such, one would except relatively larger radiative heating of the entry 
body under nonequilibrium conditions. 

Temperature distributions across the shock layer are illustrated in 
figures 11 to 13 for different conditions. For the case with no radiation, 



the heavy particle and electron temperature variations across the shock 
layer are illustrated (along with the equilibrium temperature distribution) 
in figure 11 for different entry conditions. The results show that, in the 
absence of radiation, the nonequilibrium temperature is higher than the 
equilibrium temperature throughout the shock layer for each entry condition. 
It is also noted th*t the electron temperature, which is lower than the 
heavy particle temperature during early stages of ionization, asymptotically 
approaches the heavy particle temperature during the later stages of 
ionization. As discussed in the previous section, the temperature distri- 
bution in the shock layer is relatively higher for higher altitudes because 
of higher entry velocities. For entry conditions at Z » 116 km, the 
electron temperature distributions (without and with radiation) are shown 
in figures 12 and 13 for 3 different body nose radii. As noted earlier, 
the thickness of the nonequilibrium layer decreases with increasing nose 
radius. Also, for a given nose radius, inclusion of radiation increases 
the thickness of the nonequilibrium layer. This is because the loss of 
radiation from the shock layer results in an entirely different temperature 
distribution (see fig. 13) and leaves relatively less energy for dissociation 
and ionization of the gas. 

For entry conditions at Z » 116 km, figure 14 shows the mass fraction 
of atomic hydrogen and hydrogen ion along the stagnation streamline in the 
precursor region. While equilibrium results indicate that only 5 percent 
hydrogen is dissociated and 0.018 percent is ionized, the nonequilibrium 
results show that 15 percent hydrogen is dissociated and 0.8 percent ionized. 
It should be pointed out that the composition of the precursor gas will be 
different for different entry conditions. It should be emphasized here 
again that in investigating the precursor region flow properties and their 
influence on the shock-layer flow phenomena, the entire precursor- shock 
layer solutions are obtained by iterative procedures. 

For the case with radiation and for entry conditions at Z * 116 km, the 
heavy particle and electron temperature variations across the shock layer 
are illustrated in figure 15 along with the equilibrium temperature distri- 
bution. In comparison with results of figure 11, it is seen that in the 
present case the nonequilibrium temperature is lower than the equilibrium 
temperature in certain portions of the shock. This is a direct consequence 
of the radiation cooling (i.e., radiation loss to the free stream) of the 
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shock layer. Also, in this case the nonequilibrium temperature is slightly 
higher than the equilibrium temperature in the vicinity of the wall. This 
is because cold gases near the wall absorb radiation from the high- temperature 
region of the shock layer. As would be expected, precursor heating results 
in a slightly higher shock-layer temperature distribution. 

Variations of temperature, pressure, and density along the stagnation 
streamline in the entire shock layer-precursor zone are illustrated in figures 
16 to 19 for different conditions. These results show that precursor effects 
are higher for the nonequilibrium conditions. This, however, would be expected 
since, in this case, the radiative heat flux toward the precursor region is 
considerably higher. The shock-layer nonequilibrium condition significantly 
influences the temperature and pressure variations in the precursor zone, but 
its effects on density changes are quite small. As noted earlier, in the 
shock layer, rcnequi librium results approach the corresponding equilibrium 
values at about 25 percent of the shock- layer thickness from the shock wave. 

For the equilibrium case, the influence of precursor heating on shock- layer 
temperature, pressure, and density variations is discussed in reference 2. 

For a comparison of the shock- layer flow phenomena for the two nominal 
compositions of the Jovian atmosphere, illustrative results were obtained for 
entry conditions at Z * 116 km. Results tor the temperature variation 
immediately behind the shock and for the radiative heat flux across the 
shock layer are illustrated in figures 20 to 22. It is evident from figure 
20 that the shock temperature is lower by about 2 percent for case of the 
89 percent hydrogen atmosphere. This is because, in this case, relatively 
more energy is required to dissociate the molecular hydrogen. Since the 
shock temperature is lower in this case, the radiative heat fluxes 
(q* as well as q~) are lower for both equilibrium and nonequilibrium 
conditions (see figs. 21 and 22). 

To investigate the extent of heating on an entry body, the variations 

in radiative heat flux in the shock layer were calculated for different 

conditions. As discussed earlier, the chemical nonequilibrium effects are 

more important with small body nose radius and for higher altitude entry 

conditions. Results for radiative flux toward the shock and body are shown 

in figure 23 for R* * 12 cm and 2 * 116 km. The results indicate that, in 
n 
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the nonequi librium case, the radiative heat flux is increased to about 70 
percent toward the body and almost 2.5 times toward the shock (i.e., toward 
the precursor region) . Results for radiative heating of the body for 
R* « 23 cm and Z * 143 km are shown in figure 24. The results show that the 
heat flux is about three times higher for nonequilibrium conditions. This 
is a direct consequence of the higher temperature in the nonequilibrium 
layer near the shock. 

To investigate the influence of precursor heating on viscous, nonequi- 
librium, shock-layer flow phenomena, specific results were obtained for the 
peak heating entry conditions and for an entry body with a nose radius of 
R* * 23 cm. These are presented here as final results of the present study. 

The radiative heat flux from the shock layer toward the shock front and 
the precursor region is shown in figure 25 for both equilibrium and nonequi- 
librium conditions. The results clearly indicate that heat flux toward the 
precursor region is considerably higher for nonequilibrium conditions. This 
is again a direct consequence of higher nonequilibrium temperature in the 
shock layer. As discussed before, precursor heating results in a higher 
radiative flux at the shock front. The results of figure 25 indicate that 
precursor heating results in a 15 percent increase in radiative flux in the 
nonequilibrium case while only 8.5 percent increase is noticed for the 
equilibrium condition. 

The results of equilibrium and nonequilibrium radiative flux toward the 
body (along the stagnation line) are illustrated in figure 26. Although it 
is realistic to calculate the radiative flux based on the electron temperature, 
results (for the case with no precursor effects) have been obtained also by 
ising the heavy particle temperature only for comparative purposes. The 
nonetui librium results are seen to be significantly higher than the equilibrium 
results. This is primarily due to the high- temperature region near the shock 
where nonequilibrium temperature overshoots occur. 

Figure 27 shows the variation of radiative and convective flux with 
distance along the body surface. The radiative as well as convective heat 
transfer to the body surface is seen to be enhanced by the nonequilibrium 
conditions. As discussed above, the increase in radiative heating is a 
direct consequence of higher electronic temperature. For the case with no 
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radiation, the convective heat flux toward the body was found to be the same 
for equilibrium and nonequilibrium conditions. For the case with radiation, 
however, figure 27 shows that the convective heat flux for the nonequilibrium 
case is about 20 percent higher than the corresponding equilibrium value at 
the stagnation point. This is because the cold gas near the wall absorbs 
higher radiative flux from the shock layer under the nonequilibrium condi- 
tions. As discussed before, the influence of precursor heating is enhanced 
due to nonequilibrium conditions. Figure 27 shows that precursor heating 
results in a 10.5 percent increase in the radiative flux at the stagnation 
point in the nonequilibrium case while only about 7 percent increase is 
noted for the equilibrium case. 

For the entry conditions considered in this study, therefore, it is 
logical to conclude that nonequilibrium heating of the body is significantly 
higher than equilibrium heating. Results similar to this were also obtained 
by Grose and Nealy (ref. 10) for Venusian entry conditions. For certain 
Jovian entry conditions, results presented in references 3 to 6 indicate 
that nonequilibrium heating is considerably less than the equilibrium 
heating. This obviously is in contradiction to the present findings. It 
should be pointed out that for the entry conditions considered in this study, 
the temperature just behind the shock is very high and all hydrogen molecules 
are completely dissociated. Under these conditions, H fl line emissions are 

P 

higher than for the equilibrium conditions. This is because, in addition to 
high temperature, the number density of atomic hydrogen is considerably 
higher than the equilibrium value. Thus, findings of the present study 
appear to be completely justified. 

5. CONCLUSIONS 

The results of present investigation reveal that there exists a nonequi- 
librium layer of considerable thickness in the shock layer. Since dissocia- 
tion reactions occur only within a very short distance from the shock wave, 
it is reasonable to assume that all hydrogen molecules are dissociated (half 
frozen condition) immediately behind the shock. The thickness of the nonequi- 
librium layer increases with increasing altitudes, decreasing body ncse radius, 
and with inclusion of the radiative heat-flux term in the energy equation. 
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Under nonequilibrium conditions, temperature (heavy particle as well as 
electronic) overshoots occur near the shock wave. As a result of this, the 
radiative as well as convective heat transfer to the body surface is increased 
significantly. The influence of precursor heating is enhanced due to nonequi- 
librium conditions: a 9.5 percent increase in the stagnation point radiative 
heating has been observed at an altitude of 116 km. 
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Figure 2. Species concentration variation in nonequilibrium 
shock-layer region for £ » 0. 




V * y/y $ * y*/y i 


Figure 3. Equilibrium and nonequilibrium v-velocity component and density 
variation in the shock layer for Z 3 0. 
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Figure 5. Shock-standoff variation with distance along 
the body surface. 


25 




Figure 6. Species concentration in the 
£ * 0 (with no radiation) . 
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Figure 7. Species concentration in the shock layer for 
C * 0 (with radiation). 
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Figure 9. Species concentration for different body radius in the shock 
layer for 5=0 at 2 « 116 km (with no radiation). 
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Figure 10(a). Variation of absorption coefficient for equilibrium 
and nonequilibrium conditions in the shock laver, 

Z * 116 km. 
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Figure 10(c). Variation of absorption coefficient for equilibrium 
and nonequilibrium conditions in the shock layer. 

2 » 116 km. 






Figure 11. Heavy particle and electron temperature variation 
in the shock- layer nonequilibrium region. 
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Figure 12 


. Temperature variation for different body nose 
radius in the nonequilibrium region at C « 0 
(with no radiation) . 








Figure 13. Temperature variation for different body nose radius in 
the nonequilibrium region at £ ■ 0 (with radiation). 
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Figure 14. Species concentration in the precursor region 
for 5*0. 
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Figure 15. Equilibrium and nonequilibi’um temperature variation 
in the shock layer (with radiation). 
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Figure 17. Temperature variation in the shock/precursor region along 
the stagnation streamline, Z * 143 km. 
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Figure 18. Pressure variation in the shock/precursor region 
along the stagnation streamline. 
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Figure 19. Density variation in the shock /precursor region along che stagnation streamline. 
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Figure 20. Equilibrium and nonequilibrium shock temperature variation 
for different free gas stream compositions. 
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Figure 21. Equilibrium radiative heat flux variation in the shock 
layer for different free-stream gas compositions. 
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Figure 22. Nonequiiibrium radiative heat flux variation in the shock 
layer for different free- stream gas compor ’ tions . 
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Figure 25. Equilibrium and nonequilibrium radiative heat flux towards the shock for 
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Figure 27. Variation of radiative and convective heat flux with distance along *the body surface 













